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Abstract 

The propagation of zero sound in a spin-polarized Fermi gas under harmonic con- 
finement is studied as a function of the mean-field interactions with a second Fermi 
gas. A local-density treatment is compared with the numerical solution of the 
Vlasov-Landau equations for the propagation of density distortions in a trapped 
two-component Fermi gas at temperature T = 0.2 Tp. The response of the gas to 
the sudden creation of a sharp hole at its centre is also studied numerically. 
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1 Introduction 



Current experiments on the collisional properties of a Fermi gas in two spin- 
polarized states under confinement have proven that this setup provides an 
important tool for investigating the dynamics of dilute quantum gases. For 
such a system in harmonic trapping the existence of zero and first sound 
excitations has been demonstrated experimentally [1,2] and numerically [3] by 
measuring the damping of the oscillations of a non-equilibrium density profile 
as a function of the collision rate. These studies have followed the evolution of 
the dynamical behaviour of the quantum gas from a collisional (first sound) 
regime to a collisionless (zero sound) regime. However, no direct measurement 
of the speed of sound propagation has been carried out so far in confined Fermi 
gases, contrary to the case of Bose-Einstein condensates in elongated traps [4] 
where the propagation of density distortions has been followed experimentally. 
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The main purpose of the present study is to evaluate by both theoretical 
and numerical means the propagation of density distortions in the zero-sound 
regime for a degenerate Fermi gas contained inside a cigar-shaped harmonic 
trap, in situations that are accessible to experiment. It may be useful to recall 
at this point the earlier studies of these collective excitations in homogeneous 
Fermi fluids. The experiments of Abel et al. [5] on acoustic waves propagating 
through liquid 3 He demonstrated that the transition from collisional to colli- 
sionless dynamics is signalled by a drastic change in the frequency dependence 
of sound-wave attenuation and showed that the difference in the speeds of first 
and zero sound is quite small. This is as expected for a strongly coupled Fermi 
liquid from the Landau theory [6] . The evolution of the zero-sound dispersion 
relation in liquid 3 He into the microscopic regime was determined by Skold et 
al. [7] in inelastic neutron scattering experiments and accounted for theoreti- 
cally by Aldrich et al. [8] within a polarization potential approach extending 
Landau's Fermi-liquid theory to larger momenta and energies. Landau's the- 
ory has also been used by Yip and Ho [9] to evaluate the collisionless excitation 
modes corresponding to density fluctuations and to coherent interconversion 
of various spin species in a homogeneous Fermi gas. 

At variance from the strongly coupled 3 He liquid, in a dilute Fermi gas the 
interactions can be expected to introduce only slight modifications to the 
bare sound velocities, which in the homogeneous noninteracting limit are vf 
and vf/V3 for zero and first sound, respectively. Here, the Fermi velocity vp 
is determined by the Fermi wave number and hence by the particle number 
density. These well-known results [6] immediately suggest how the propagation 
of collective density fluctuations in a confined Fermi gas may be treated within 
a local-density approximation. 

In this work we first present an analytical study of the behaviour of the 
Fermi velocity and of the speed of zero sound in a two-component inhomoge- 
neous Fermi gas as functions of the mean-field interactions in the collisionless 
regime. We support our analytical estimate by numerical simulations solving 
the Vlasov-Landau equations (VLE) for the fermionic Wigner distributions at 
T = 0.2 Tp, Tp being the Fermi temperature [3]. We evaluate the propagation 
of zero sound by simulating an experiment where perturbation in the density 
profile is introduced by a laser beam permanently focused on the centre of the 
trap [10,11]. We also simulate an experiment carried on by Dutton et al. [12] 
on a Bose-Einstein condensate, where a defect with a large gradient in the 
density is generated by cutting off a slice of the profile by means of an ultra- 
compressed slow pulse of light. In their experiment the generation of quantum 
shock waves in the condensate was observed. 

The paper is organized as follows. In Sec. 2 we recall the equations of motion 
for a two-component Fermi gas and evaluate in a semiclassical local-density 
approach the Fermi velocity and the zero sound velocity for such a system. In 
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Sec. 3 we report illustrations of our numerical simulations for a one-component 
and a two-component Fermi gas, and compare the results with our theoretical 
estimates. Conclusions and perspectives are given in Sec. 4. 



2 The model 



The two spin-polarized fermionic components of the g described by the 
Wigner distribution functions f^(r,p,t), with j — 1 or 2. These obey the 
kinetic equations 

^ + — • V/&'> - VC/W • Vp/&'> = Ci 2 [/W]. (1) 

At high dilution the interactions between atoms of the same spin are sup- 
pressed by the Pauli principle, so that the mean-field Hartree-Fock (HF) ef- 
fective potentials are U^\r,t) = V^.}(r,t) + gn®(r, t) with j indicating the 
species different from j and v£j}(r,t) being the sum of external confining and 
driving potentials. Here we have set H — 1 and g = 2na/m r , with a the s-wave 
scattering length (a > for the present case of repulsive interactions) and m r 
the reduced mass. The total spatial densities (r, t) = (r)+8n^ (r, t) are 
the sum of the equilibrium densities (r) and the time- dependent perturba- 
tions 8n^\r,t), and are obtained by integration of f^ j \r,p,t) in momentum 
space. Finally C\2[f^} is the collision integral, which will be set to zero in 
the collisionless regime. Equation (1) in this regime will be solved numerically 
in Sec. 3, using the algorithm developed by Toschi et al. [3] and taking the 
equilibrium density profiles as given in the HF approach [13] by 



ni j) ( 



r >=/^MKfi +aW(r) -4 



(2) 



where (3 = 1/ksT and lft> is the chemical potential of species j. 



In the rest of this section we present instead a theoretical discussion of the 
long- wavelength excitations of the gas. The perturbation Sn^ (r, t) moves with 
a group velocity whose natural unit is the Fermi velocity Vp\ In an inhomo- 
geneneous gas the Fermi velocities depend on position and will be related 
below to the local chemical potential of each species through a local density 
approximation. Assuming for the moment that they are known, the dynamical 
response of the two-component Fermi gas can be calculated within a random 
phase approximation [14,15] from the set of eigenvalue equations 

8n® (q, fl) = - Xj (q, A) [5V® (q, Q) + g 5n® (q, Q)] , (3) 
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where 



with M 



Xj = Mj 



2v { pq 



In 
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n-v<i ) q / 



(4) 
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m 2 /27r 2 being the density of states of each component per 
unit volume at the Fermi level. In writing these equations we have specifically 
assumed that we are treating the case of long-wavelength excitations propa- 
gating with wave number q along the axis of an elongated cylindrical trap, 
with the Fourier transform 8V^\q, Q) of the driving potentials being inde- 
pendent of the radial coordinate and the Fourier transform Sn^ (q, f2) of the 
density distortions being taken at the value r = of the radial coordinate. 
Following [9] and [15] it is easily seen that the collective modes of the Fermi 
gas described by Eq. (3) are determined by the solutions of the equation 



i-g 2 xi(q,n) X 2(q,n) = o 



(5) 



yielding in general two eigenfrequencies £lj(q). 



Let us turn next to the calculation of the Fermi velocities. Each of them is 
determined by the corresponding single-particle dispersion relation as vp = 
dE(k)/dk\k=k F , yielding in a local density approximation 



Ji) 
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where /4oc( r ) = ^ ~~ U^\r) is the local chemical potential. In cylindrical 
harmonic confinement and in the absence of interactions v$ can be written 

as 



(r, z) = [2 (ji® - V {j) (r, z)) frrij] 1/2 (7) 

where V^(r,z) = rrijtUj(r 2 + e'jZ 2 )/2, with radial trap frequency ujj and 

anisotropy Sj. Under driving the radial and axial components of vp^ then 
obey harmonic equations of motion, 



j { p(t)\ z=0 = (2^'V^') 1/2 sin(^t + r ) 



(8) 



and 



vf{t)\ r= ,= {2^/m, 



1/2 



sm(e jujjt + 2 ), 



(9) 

where the phases <p r and <p z are determined by the initial value of the pertur- 
bation. In the presence of mean-field interactions, on the other hand, Eq. (7) 
is replaced by 



r,z) = 



U) _ v®(r, z) - gn£\r, z) ) /m 



1/2 



(10) 



In this case harmonicity under driving will be lost, except when the mean field 
term can be treated perturbatively and for small axial displacements from to 
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the center of the trap. In this regime it is easy to demonstrate that the mean- 
field interactions lower the Fermi velocity. More generally, their net effect is to 
decrease the chemical potential by inducing a decrease of the fermion density. 

With the above as background, we conclude this section by presenting our 
estimates for the speed Cq of collisionless sound waves propagating along the 
trap axis at r = in two special cases, that will be tested in the next section by 
comparison with the evolution of density profiles in time from the numerical 
solution of Eq. (1): 

(i) A single- component Fermi gas. Setting g — 0, Eqs. (3) and (4) immediately 
yield 

c o = v F (0,z). (11) 

This situation will arise not only in the obvious case where the interspecies 
mean-field interactions vanish, but also in the case where they become so 
strong in a two-component Fermi gas that the regime of phase separation pre- 
dicted by the theory [13] is being approached. In this latter case the residual 
mean-field interactions with the second species merely shift the chemical po- 
tential of the species that resides in the inner part of the trap and hence will 
shift its Fermi velocity. 

(ii) A symmetric two- component Fermi gas. Taking equal masses and numbers 
of particles for the two fermionic species inside identical traps, we can write 
c = Xvp for the local speed of sound wave propagation and determine A 
from Eq. (5) using the Landau Fermi-liquid theory to describe the mean-field 
interactions through the appropriate Landau parameters. Only the spin up- 
down scattering is allowed in the s-wave channel for a two-component mixture 
of spin polarized fermions, so that the only non-zero interaction parameters 
are Fq and F a with the condition Fq = — F a (= F , say). Equation (5) can 
be written in the form 

-TO 

with 

F = gj\f = 2akp/n. (13) 

In Figure 1 we have plotted the behaviour of the solution of Eq. (12) as a 
function of the Fermi velocity for various values of the scattering length. We 
also show in Figure 2 the values of vp and c as functions of the scattering 
length, for the choice of the other system parameters that will be adopted 
in the numerical calculations reported in the next section. It is seen from 
these figures that, in order to obtain a very significant increase of the speed 
of zero sound over the Fermi velocity from repulsive mean-field interactions, 
one would need to work on a high-density system having huge values of the 
scattering length. The value of a can to some extent be tuned by exploiting 
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Feshbach resonances [16], but strong repulsions will also favour demixing of 
the two-component Fermi gas. 



3 Numerical results 

3. 1 Perturbation introduced by a permanent laser beam 

As a first case we consider a one-component gas of 1000 spin-polarized 40 K 
atoms inside an axially symmetric trap with u = 100 s _1 and e = 0.1. We 
use 200 computational particles to represent each fermion, in order to obtain 
a low noise-to-signal ratio [3]. The equilibrium density profile is perturbed 
by an effective potential U exp(— z 2 /2a 2 ) with amplitude Uq and width a, 
which simulates a laser beam acting continuously on the centre of the trap. 
The fermions interacting with the laser are pushed out of the trap centre 
and the density distortion moves towards the ends of the trap with a ve- 
locity which corresponds to the local Fermi velocity given by Eq. (7). Fig- 
ure 3 shows the evolution in time of the perturbation along the z axis for 
r = 0. For such a system the maximum value of the axial zero-sound velocity, 
which is equal to the axial Fermi velocity, can be evaluated using Eq. (7) as 
c^ aa; | r= o = v^ ax \ r= o = 0.383 aho/ms, where a^o = (h/muj) 1 / 2 is the radial har- 
monic oscillator length. Figure 3 also shows that the positions of the density 
peaks as calculated numerically at various time steps are in good agreement 
with those estimated from Eq. (7). 

We consider next a 50-50 mixture of 2000 40 K atoms which are polarized in the 
7/2 and 9/2 spin states, under the same external potentials as in the previous 
case. The value of the scattering length of 40 K (a = 80 Bohr radii a ) is so small 
that vf is lowered by the interactions by only 0.1%. From Fig. 2 it is evident 
that the speed of zero sound is practically unchanged. To see the interaction 
effects we study the case a = 2 x 10 4 ao- For this value of the scattering 
length the Fermi velocity is obtained from the numerical calculation of the 
chemical potential and of the effective potential as i^ o:E | r=0 = 0.343 aho/rns. 
The corresponding values of A and Co are 1.01 and 0.347 a^/ms respectively 
(see Fig. 2). The time evolution of the perturbation in the two components of 
the interacting Fermi gas is shown in Fig. 4, in comparison with that of the 
noninteracting gas. 

At still higher couplings spatial phase separation starts in the two-component 
Fermi gas. The case a = 8 x 10 4 a is illustrated in Figure 5. The two compo- 
nents are amost completely demixed and only one of them feels the external 
perturbation. In this case the interactions with the second component merely 
shift the chemical potential entering the evaluation of the Fermi velocity of 
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the first component, which equals the speed of propagation of zero sound. 
3.2 Perturbation introduced by a compressed slow light pulse 

In the work of Dutton et al. [12] ultra-compressed slow light-pulses were used 
to pump a small slice of a Bose-Einstein condensate into an untrapped state, 
leaving behind a sharp hole in the density profile. The density defect in the 
condensate evolves through the formation of solitons, this process being anal- 
ogous to the creation of shock waves in a classical fluid. 

We generate the same type of initial profile in a two-component Fermi gas for 
various shapes of the density vacancy (V-shape or U-shape of various sizes) 
and examine the effects of the density gradient and of the size of the defect. No 
coherent behaviour is observed in the time evolution of the density profiles. In 
Fig. 6 we plot the time evolution of one of the two components of the Fermi gas 
for various values of the scattering length (a = 0, 800 a and 2 x 10 4 a ) in the 
case of a small U-shape vacancy. In the two former cases the particles in filling 
the vacancy generate some density waves, which reach a maximum amplitude 
after about 150 ms. This is the time lag nedeed for the perturbation to arrive 
at the ends of the trap and to be reflected there, as can be roughly estimated 
by dividing the axial half-length L z of the trap by c™ ax giving L z /c™ ax ~ 
140 ms. Because of the absence of thermal dissipation in our model, the density 
waves continue to travel along the profile and are again reflected at the ends 
of the trap with a periodicity of about 2L z /c™ ax ~ 280 ms. The revivals of 
smaller amplitude which are seen to occur at about 300 ms correspond to the 
situation in which the perturbation is passing through the centre of the trap. 
On increasing the coupling strength up to a = 2 x 10 4 a (right panel in Fig. 6), 
we find that the restoring force due to the other fermionic component plays a 
role in depressing the intensity of the density waves relative to those seen in 
the weak-coupling regime. 



4 Conclusions 

In summary, we have studied the propagation of collisionless sound waves in a 
trapped two-component Fermi gas, assuming a mean-field model for the repul- 
sive interactions between the two fermionic species and neglecting dissipative 
collisions as described by a Boltzmann quantum collision integral. The inclu- 
sion of such collisions is needed to examine the transition from the zero-sound 
to the first-sound regime [3] , that we expect will be characterized by a marked 
change in the speed of sound propagation. We also expect that zero-sound 
waves will be Landau-damped in the case of attractive interactions, as is the 
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case for a homogeneous Fermi liquid [6] . 

We have found that, both at weak coupling and in the strong-coupling regime 
where the two Fermi clouds are demixing from their mutual repulsions, the 
speed of zero-sound waves is very close to the Fermi velocity as determined by 
the local particle density through the local chemical potential. These results 
have been obtained from comparing semi-analytical estimates with numerical 
studies of the evolution of peaks in the density profiles under driving from an 
external potential simulating a blue-detuned laser that keeps expelling atoms 
from the central region of the trap. 

We have also examined numerically the time evolution of impulsively created 
density defects which are characterized by sharp density gradients in a two- 
component Fermi gas with repulsive interactions. In the collisionless regime we 
have found that the filling of the vacancy at the trap centre generates density 
waves which travel to the ends of the trap and are back-reflected there, so that 
they appear as periodic revivals of distortions in the density profile. These re- 
vivals become attenuated as the strength of the coupling is increased, presum- 
ably as a consequence of the increasing efficiency of momentum-redistribution 
processes through non-dissipative collisions at strong coupling. 
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Fig. 1. The A factor as a function of vf (in units of flfe/ms with a^ = (h/mtu) 1 / 2 ) 
for various values of the scattering length: a = 2 x 10 4 ao (continuous line), 8 x 10 3 ao 
(short-dashed line), 2 x 10 3 ao (dotted line) and 800 ao (long-dashed line), ao being 
the Bohr radius. 
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Fig. 2. The speed of zero sound (continuous line) and the Fermi velocity (dashed 
line), in units of a/, /ms, as functions of the scattering length a (in units of cto) for a 
two-component Fermi gas. The calculation refers to 2000 particles in an anisotropic 
harmonic trap with u> = 100 s _1 and e = 0.1 and the values of the velocities refer 
to the case of a little perturbation near the centre of the trap. 
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Fig. 3. Time evolution of a defect introduced in a one-component Fermi gas, plotted 
as a function of the axial distance z (in units of aho)- The bottom profile refers to 
the initial state and the subsequent profiles are separated by time intervals of 10 ms. 
The dots indicate the analytical estimate of the position of the perturbation peak 
from Eq. (7). 
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Fig. 4. Time evolution of a defect introduced in both components of an interacting 
two-component Fermi gas with a = 2 x 10 4 ao (dashed lines), in comparison with the 
one-component case (continuous lines). From bottom to top each profile corresponds 
to a time lag of 10 ms. Dots indicate the analytical estimate of the position of the 
peak in the non-interacting case from Eq. (7) and crosses indicate the peak positions 
corresponding to the zero sound velocity evaluated at t = from Eqs. (10), (12) 
and (13). 
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Fig. 5. Time evolution of a defect introduced in one component of an interacting 
two-component Fermi gas with a = 8x 10 4 ao in a regime of phase separation (dashed 
lines). From bottom to top each profile corresponds to a time lag of 5 ms. The dots 
indicate the analytical estimate of the position of the peak of the perturbation from 
Eq. (7). 
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Fig. 6. Time evolution of the axial density profile as a function of z (in units of a^ G ) 
for the case of a small U-shape vacancy with a = 0, 800 ao and 2 x 10 4 ao (from left 
to right). From bottom to top, the time lag of each profile (in ms) is indicated in 
the figure. The profiles which are not shown in the ranges of time lag between 200 
and 300 ms and between 325 and 450 ms are essentially smooth. 
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